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ABSTRACT 

We build a large sample of Seyfert 2 galaxies (Sy2s) with both optical spec- 
tropolarimetric and X-ray data available, in which 29 Sy2s with the detection of 
polarized broad emission line (PBL) and 25 without. We find that for luminous 
Sy2s with L[o m] > 10^^ erg s~^, sources with PBL have smaller X-ray absorption 
column density comparing with those without PBL (at 92.3% confidence level): 
most of the Sy2s with Nh <W'^^-^ cm'^ show PBL (86%, 12 out 14), while the 
fraction is much smaller for sources with heavier obscuration (54%, 15 out 28). 
The confidence level of the difference in absorption bounces up to 99.1% while 
using the "T" ratio (F2_io kev/F[o iii]) as an indicator. We rule out observation 
or selection bias as the origin for the difference. Our results, for the first time 
with high statistical confidence, show that, in additional to the nuclei activity, 
the nuclear obscuration also plays an important role in the visibility of PBL in 
Sy2s. These results can be interpreted in the framework of the unified model. We 
can reach these results in the unified model if: a) the absorption column density 
is higher at large inclinations and b) the scattering region is obscured at large 
inclinations. 

Subject headings: galaxies: active — galaxies: Seyfert — polarization 



1. Introduction 

The AGN unification model proposes that Seyfert 1 and 2 galaxies (Syls and Sy2s 
hereafter) are intrinsically the same objects, and the absence of broad emission lines in 
Sy2s is ascribed to the obscuration along the line of sight by a pc-scale dusty torus (see 
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the review by Antonucci 1993). The most convincing evidence is the detection of polarized 

broad emission Unes (hereafter PBL) in some Seyfert 2 galaxies (Antonucci & Miller 1985; 
Tran 1995; Heisler, Lumsden & Bailey 1997; Moran et al. 2000; Lumsden et al. 2001; 
Tran 2001). Similarly infrared (IR) observations showed the existence of obscured broad 
line regions (BLRs) in Sy2s (Veilleux, Goodrich & Hill 1997). Further evidence supporting 
the unification model comes from X-ray observations of Sy2s which show large amounts of 
obscuration, typically above 10^^ cm~^ (e.g., Turner et al. 1997; Bassani et al. 1999). 

Despite observations do generally support orientation based unification model, only ~ 
50% of Sy2s show broad lines in the polarized spectrum (e.g. Tran 2001, Gu & Huang 
2002). With an optical spectropolarimetric study of a well-defined and statistically complete 
IRAS 60-/xm selected Sy2 sample, Heisler et al. (1997) found a relationship between the 
detectability of polarized broad Ha and the IRAS /60//25 flux ratio that only those galax- 
ies with warm IRAS colors {fm/f 25 < 4.0) show PBL. Heisler et al. suggested that the 
detectability of PBL simply depends on the inclination of the torus to the line of sight: in 
a Sy2 with the torus highly inclined, cooler infrared color is expected, and the broad-line 
scattering screen could also be obscured. 

A simple prediction of the inchnation-related model is that Sy2s without PBL (hereafter 
NPBL Sy2s) should show higher absorption column density since they are more inclined 
than the Sy2s with PBL (hereafter PBL Sy2s). However, following studies have claimed 
no difference in the absorption column density between two types of Sy2s (Alexander 2001; 
Tran 2001, 2003). Furthermore, as Alexander (2001, also see Lumsden et al 2001; Tran 2001, 
2003; Gu & Huang 2002) pointed out, the difference in the IRAS /60//25 flux ratio is not an 
good indicator of the inclination but the relative strength of galactic and Seyfert emission. 
These studies (also see Cheng et al. 2002; Lumsden & Alexander 2001), instead, have shown 
that the presence of PBL in Sy2s depends on the AGN luminosity: Sy2s with PBL have 
higher luminosity comparing with Sy2s without PBL. Explanations to the observational 
results include: a) The contribution from the host galaxy or from a circumnuclear starburst 
would dilute the nuclear optical spectrum, making the detection of PBL more difficult for 
Sy2s with lower luminosity (Alexander 2001; Gu et al. 2001); b) Alternatively, Tran (2001, 
2003, also see Yu & Hwang 2005) suggested that at least some of the Sy2s without PBL 
are powered by starburst rather than accretion onto a supermassive black hole, therefore, 
the BLRs simply do not exist; c) More luminous sources tend to have large scaleheight of 
the scattering region thus increasing the visibility of PBL (Lumsden & Alexander 2001); 
d) Nicastro, Martocchia & Matt (2003) have argued that at very low accretion rates (and 
therefore lower luminosities) the clouds of the BLRs would cease to exist and the absence 
of PBL in Sy2s is consistent to their low accretion rates, e) In the case of low luminosity 
nuclei, the adjacent bright sources can easily outshine the nuclear flux and the Nh derived 
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from X-ray spectrum may be underestimated (e.g. Georgantopoulos & Zezas 2003). f) The 

large-scale dusty environment (Panessa & Bassani 2002) or complex and variable obscuring 
material (Matt 2000a; Risaliti 2002) may to some extent affect the appearance of PBL in 
Sy2s. g) Long term large amplitude variation in the nuclei activity could vary the PBL flux 
and thus the detectability of PBL (Lumsden et al. 2004). 

Meanwhile, it's worth to note that there are also weak evidences showing different 
absorption in two types of Sy2s: Gu, Maiolino & Dultzin-Hacyan (2001) found slightly (but 
not statistically conclusive) lower Nh for PBL Sy2s, and Lumsden, Alexander & Hough 
(2004) found a considerably higher detection rate of scattered broad in a small sample 
of Compton-thin Sy2s. Note one must be cautious while comparing the fraction of PBL 
Sy2s between samples since the luminosity might have played a major role. These evidences 
suggest that besides the AGN luminosity, the X-ray absorption column density also plays 
a role in the visibility of PBL in Sy2s. We point out that the role of the absorption (if 
exists) could reveal itself in a luminous Sy2 sample where the influence of luminosity on the 
visibility of PBL is weak enough. 

In this paper we revisit the issue of whether the nuclear obscuration in Seyfert 2 galaxies 
affects the visibihty of PBL by focusing on luminous Sy2s. The launch of the Chandra X-ray 
Observatory in 1999 and IQsilsi- N ewton in 2000 open a new era of X-ray astronomy. New 
Chandra and IQs/lls/l-N ewton observations have significantly enlarged the sample of Sy2s with 
both spcctropolarimetric and X-ray observations available, and also provide more reliable 
X-ray measurements thanks to their much higher spatial resolution and better sensitivity. 
In this paper, we present a large sample of Sy2s, for which both the spcctropolarimetric 
and X-ray data are available, to probe the nuclear obscuration for PBL/NPBL Sy2s. Our 
sample consists of 29 PBL Sy2s and 25 NPBL Sy2s. Among them 8 Chandra and 30 XMM- 
Newton observations are available either from literature or from archive. Throughout this 
paper we use the cosmological parameters Hq = 70 km s~^ Mpc~^, = 0.27 and Q,\ — 0.73. 



2. Sample selection 

We collect all Sy2s with spcctropolarimetric observations from literature (from 1985 
to 2006, see table 1). The spcctropolarimetric data are mainly from several large surveys 
including: the infrared-selected sample of Hcislcr ct al. (1997), the far infrared flux and 
luminosity limited sample of Lumsden et al. (2001), the distance limited sample of Moran 
et al. (2000, 2001), the heterogeneous optical and mid infrared selected sample of Tran et 
al. (2001), and infrared color selected sample of Young et al. (1996). We then exclude NGC 
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2992, NGC 5506, NGC 5252, NGC 7314, MCG -3-34-64 and Mrk 33^1 from discussion due 
to their intermediate classification (i.e. Syl.Ss, 1.9s) in NED. To avoid luminosity selection 
bias due to the redshift difference when comparing properties between the PBL and NPBL 
Sy2s, we confine our sample within z < 0.06. We collect [O ill] A5007 and X-ray data 
from literature, and present spectra analysis of archive Chandr a /XMM- Newton data for 
8 sources in §3. The result from most recent observation is adopted when two or more 
observations exist. Note in the table, there are 8 upper limits in the hard X-ray flux due 
to X-ray non-detection in the hard band and their X-ray absorption column densities in 
literature were estimated either from soft band X-ray data or the strength of their X-ray 
emission relative to optical band. We exclude these Nh from our following analysis. NGC 
4117 is also excluded since its [O ill] A5007 flux is not available from literature. This leaves 
a sample composing of 2s[^ Sy2s with PBL, and 25 Sy2s without PBL, for which both the 
spectropolarimetric and X-ray data are available. Fig. 1 shows the redshift against the 
luminosity of extinction-corrected [O ill] A5007 emission for PBL and NPBL Sy2s in the 
sample. As previous studies have shown, we clearly see higher luminosities for PBL Sy2s 
(with a confidence level of 99.7%), indicating the nulear activity plays a major role in the 
visibility of PBL in Sy2s. 

The optical and X-ray data are presented in Table 1. The table lists, in turn, the 
name of the galaxy, redshift z as reported in NED, the spectropolarimetric properties, the 
corresponding references, the extinction-corrected flux of [O ill] A5007 emission in units of 
10~^^ erg s~^ cm~^, the luminosity of extinction-corrected [O ill] A5007 emission in units of 
erg s~^, the references for the [O ill] A5007 emission, the observed rest frame hard X-ray 
(2-10 keV) flux in units of erg s~^, the X-ray absorption column density (Nh) in units of 
cm~^, the equivalent width (EW) of the fluorescence iron line in units of eV, the references 
for the X-ray properties. The luminosity of extinction-corrected [O ill] A5007 emission is 
given as L[o mj = 47rD^F|^Q''jjjj, where F|^q''jjjj is the extinction-corrected flux of [O ill] A5007 
emission derived from the relation (Bassani et al. 1999) 

TT'COr 77'Obs (HQ,/H/3)obsN2.94 

We assume an intrinsic Balmer decrement (Hq,/H^)o = 3.0. 



^Among the six sources, NGC 5506 and Mrk 334 do not show PBL in the polarized spectra and the rest 
four do. Including these sources in our sample does not make significant difference to our major results. 

^For one galaxy IRAS 04385-0828, the value of Nh is unavailable, but we can get the reference to hard 
X-ray flux from PoUetta et al. (1996). When using the "T" ratio for analysis, the number for PBL Sy2s is 
30. 
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3. X-ray spectral analysis 

In this section, we report the results of X-ray spectral fitting to archive Chandra and 
'yjsHyi- Newton spectra of 8 Sy2s in the sample. The data were reduced using CIAO 3.2.2 
and XMMSAS 6.5.0 respectively. The size of each source on the detector was estimated in 
order to determine appropriate source extraction regions, typically ~ 2" radius {Chandra) 
or ~ 30" {XMM-Newton) for on-axis point sources. The background spectra were extracted 
from source-free annulus around the source. The spectra of each galaxy were binned to a 
minimum of 1 counts per bin and we adopt the C-statistic (Cash 1979) for minimization. 
Spectral fits were performed using XSPEC version 11.2 in the 0.5 - 8 keV band. All the 
quoted errors are 90% confidence range for one parameter of interest. 

Each spectrum was initially fitted with a simple model consisting of a powerlaw plus 
Galactic and intrinsic neutral absorption. In many cases this simple parameterization is not 

sufficient to model the whole 0.5 - 8 keV spectrum. Residuals often show a soft excess on 
top of the powerlaw. The soft excess is fitted here as a scattered powerlaw component (with 
the same powerlaw slope but no intrinsic absorption). The possible presence of a narrow 
emission line centered at 6.4 keV originating from neutral iron has also been checked, and 
modeled with a single Gaussian line. 

We note that in Compton-thick sources with > 10^^ cm^^, the transmitted compo- 
nent is heavily suppressed below 10 keV and the spectrum observed in the 2-10 keV band 
might be dominated by the reflection component (Matt et al. 2000b). In this paper, NGC 
34, NGC 3982, Mrk 573, Mrk 1066 are classified as Compton-thick based on their large Fe 
Ka EW (> IkeV except for NGC 34, see notes in Appendix) and small F2_io kev/F[o iii] 
ratios (<0.1) (Maiolino et al. 1998, Bassani et al. 1999, Guainazzi et al. 2005b). For 
Compton-thick sources, we use the refiection model {pexrav model in XSPEC; Magdziarz & 
Zdziarski 1995) for spectrum fitting. 

Given the purpose of this work (to obtain a proper description of the spectra in terms 
of absorption, 2-10 keV fiux, and Fe K fine intensity), these simple parameterizations yield 
adequate fits to all the spectra presented here. The best-fit spectral parameters are listed in 
Table 2 and notes on individual objects are given in Appendix. 

4. different obscuration in two types of Sy2s 

In Fig. 2 (b), we plot the luminosity of extinction-corrected [O ill] A5007 emission versus 
Nh- The separation is apparent for two types of Sy2s. The diagram can be roughly divided 
into three regions with the boundaries at Nh = 10^^-^ cm~^ and L[o iii] = 10^^ erg s~^. For 
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the luminous Sy2s (with L[o m] > lO^^erg s"^), we can clearly see that most of the Sy2s with 
Nh <10^^'^ cm~^ show PBL (86%, 12 out 14), while the fraction is much smaller for sources 
with heavier obscuration (54%, 15 out 28). For the Sy2s with lower [O ill] luminosity (< 
IQ^i erg s~^), only a small fraction show PBL (17%, 2 out of 12), and due to the limited 
number of sources, we are not able to tell if the fraction of PBL sources depends on the 
X-ray absorption at lower luminosity. 

In Fig. 3 (left panel) we plot the Nh distributions for all Sy2s with/out PBL. Since 
there are 11 censored data (lower limits) among PBL Sy2s and 10 among NPBL Sy2s, we 
use the survival analysis methods ASURV (Fcigclson & Nelson 1985) for statistical analysis. 
We find little difference (with a confidence level of 66.5% of the difference, see table 3) in Nh 
between PBL/NPBL Sy2s and the mean values of log Nh (in units of cm~^) are 23.755±0.19 
and 23.852 ± 0.274, respectively (for NGC 4501 and NGC 7590 we adopt the Nh upper 
limits as the measured values since ASURV could not deal with the case which contains 
both upper and lower limits). However, if we only consider the luminous Sy2s with L[o iii] 
> 10^^ erg s~^, K-S test shows that the probability for two samples to be extracted from the 
same parent population is about 7.7%, and the mean values of log Nh are 23.739 ±0.212 and 
24.428 ± 0.192, respectively (Fig. 3, right panel). The results suggested that for luminous 
Sy2s in our sample, sources without PBL show larger obscuration than those with PBL with 
a confidence level of 92.3%. 

To further examine if obscuration plays a role in the detection/visibility of PBL in 
Sy2s, we explore other potential measures of obscuration. By studying a large sample of 
Sy2s, Bassani et al. (1999) found that the "T" ratio F2-10 kev/F[o iii] is a good indicator of 
nuclear obscuration. In particular, it is anticorrelated with both the column density Nh and 
the Fe Ka line EW, and that these quantities can be used as probes of the obscuration to 
the center of the AGN. In Fig. 2 (c, d) we plot the Fe Ka line EW and "T" ratio versus the 
luminosity of extinct ion- corrected [O ill] A5007 emission. Similar patterns as seen in Fig. 2 
(b) are also obvious that for luminous Sy2s, the NPBL sources tend to be more obscured. 
Fig. 4 shows the "T" ratio distributions for all Sy2s (left) and for luminous objects only 
(right). A K-S test shows that the possibility for these two samples to be extracted from 
the same parent population is about 25%. The mean value of the log T are -0.087±0.145 
and -0.342±0.217 for Sy2s with and without PBL respectively. Similarly the confidence level 
for the difference is much higher (at level of 99.1%) for luminous Sy2s only. Turning to 
the Fe Ka line EW, K-S tests also confirms that for luminous Sy2s the difference between 
the two sample is present (at 95.3% level) with available data. The mean values of log 
EW(Fe) are 2.626±0.107 and 2.999±0.066 respectively. After examining three independent 
indicators for obscuration, we conclude that for luminous Sy2s, sources without PBL have 
higher obscuration than those with PBL, confirming the suggestion that the obscuration does 
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play an important role in the detectability /visibility of PBL. The results from K-S tests and 
average values for Nh, "T" ratio and FeK line EW are summarized in table 3. 

5. Discussion 

It's clear that, as many previous studies have shown, PBL Sy2s have higher luminosities 
than NPBL Sy2s (see Fig. 1 & 2), indicating the primary determinant of PBL visibility is 
the nuclei luminosity. In this paper, by focusing on luminous Sy2s, we find that the nuclei 
obscuration also plays an important role in the visibility of PBL. For Sy2s with L[o mj > 
]^q4i gj^g g-i Q^j, sample, we find that NPBL Sy2s have higher X-ray column density than 
PBL Sy2s at a significant level of 92.3%. While using the "T" ratio or the Fe K line EW 
as indicator of nuclear obscuration, the confidence level of the difference in obscuration gets 
even higher (99.1% and 95.3% respectively). Our results are consistent with Lumsden et al. 
(2004) who reported higher detection rate of PBL in Compton-thin Sy2s, but with much 
higher confidence level. Consistent with previous studies, most (83%) of the less luminous 
Sy2s (L[o III] < 10^^ erg s~^) do not show PBL, the nature of which is still unclear (see §1) and 
is beyond the scope of this paper. We also demonstrate that since most of the less luminous 
Sy2s do not show PBL independent of absorption, adding them to the sample weakens the 
difference in obscuration found in the luminous sample. This explains why previous studies, 
which did not exclude less luminous sources, found no difference in absorption. 

5.1. Selection effect? 

It is worth stressing whether the difference in the nuclear obscuration between two types 
of Sy2s could be due to possible observational and sample selection bias. We note our sam- 
ple is an amalgamation of different observations with diverse quality of spectropolarimetric 
data. We first examine whether the non-detections of PBL in the sample are due to the 
weakness/lack of PBL or due to the limited sensitivities of the spectropolarimetric data. 
Lumsden et al. (2001) showed that the S/N in their sample is sufficient for all but 2 of 
the NPBL S2s and attributed the nondetections to significantly weaker scattered fiux. Tran 
(2003) pointed out that the distributions of [O ill] fiux, which is a good indicator of the 
strength of the Seyfert nucleus, are virtually the same between two types of Sy2s in his 
sample, suggesting that the nondetections are not likely due to the detection limit of the 
survey. The sensitivity of Moran's sample (2001) is found to be even better than that of the 
other samples mentioned above (Gu & Huang 2002). In Fig. 5, we plot the "T" ratio vs. 
the extinction corrected [O ill] fiux for our composite sample. We can clearly see that for 



-8- 



luminous sources (with L[o mj > lO^-*^ erg s~^) in our sample, there is no difference in the 
[O III] flux distributions between PBL and NPBL Sy2s. For comparison, sources with lower 
luminosities are also plotted, most of which are much weaker in the [O ill] flux. We conclude 
that most of the nondetections for our luminous sources are likely due to the weakness or 
lack of PBL but not due to the limited sensitivities of the spectropolarimetric data. 

We then verify if our results are biased by combining samples with different selection 
criteria and survey depths into one single sample, i.e., if some of samples tend to select more 
obscured sources but with poor spectropolarimetric data, and/or some others tend to select 
less obscured sources but with better spectropolarimetric data. By plotting in Fig. 6 the 
"T" ratio vs. [O ill] luminosity for sources in each subsamplc, we can see this is not the case 
for our composite sample. We find that each subsample spans a similar obscuration range 
with that of the composite sample, and the difference in the obscuration between PBL and 
NPBL Sy2s are also visible in most of the subsamples. 

We made an additional test to check whether the distributions of z and [O ill] luminosity 
for our luminous sample are different (see Fig. 1 and 2 (a)). Using ASURV, we get average 
values of (z) = 0.02± 0.003 and (logL[o m]) = 42.198 ± 0.121 for PBL Sy2s whereas (^) 
= 0.017± 0.004 and (logL[o m]) = 41.991 ± 0.157 for NPBL Sy2s. The distributions of z 
and [O III] luminosity for luminous sources are similar at levels of pnuu = 29.4% and pnuu 
= 30.6%, respectively. The similarity indicates that the difference in absorption could not 
be biased by different redshift/luminosity which both affect the visibility of PBL. Also, the 
dilution effect, which is dependent of the rcdshift and luminosity, might bias the visibility of 
PBL for less luminous Sy2s, but itself alone can not explain the difference in the absorption 
between PBL and NPBL Sy2s. Actually we note that the dilution effect to the visibility of 
PBL is much weaker for our luminous sources. This can be seen from the fact that PBL 
Sy2s can be detected in most of luminous Sy2s with smaller obscuration (Nh < 10^^-^ cm~^ 
or "T" ratio > 10~°-^). The dilution effect in X-ray (to the measurement of N/^) is also 
much weaker for luminous Sy2s. We conclude that there is no observational bias which can 
produce the difference in absorption between two types of Sy2s in our sample, and a physical 
link between the visibility of PBL and nuclear obscuration is required. 

5.2. physical explanation to the difference in absorption 

The results presented here for luminous Sy2s can be interpreted within the context of 
the unified model for Seyfert galaxies, in agreement with the torus geometry portraited by 
Heisler et al. (1997): the main electron scattering is confined to a conical region that is close 
to the thickness of the torus. More inclined Sy2s could have the broad line scattering screen 
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also obscured thus make PBL weaker or nondetectable. According to the unified model, 
more incUncd sources expect heavier obscuration, thus explain the difference in obscuration 
between PBL and NPBL Sy2s. 

Note the high detection rate (86%) of PBL in luminous Sy2s with Nh < 10^^-^ cm-^ 
suggests that in most sources the Compton-thin X-ray obscuring material can not have much 
larger scale than the scattering screen, supporting the torus scheme of the unified model. 
Furthermore, while extended obscuration from the host galaxy might explain the absence of 
PBL in some sources (such as NGC 5506, see Lumsdcn ct al. 2004), it could not be the major 
cause, otherwise we should have seen a large number of NPBL Sy2s among luminous Sy2s 
with Nfj < 10^^'^ cm~^. It's interesting to note that NGC 5506 is one of two intermediate 
Seyferts without PBL detected (the other one is Mrk 334, see §2), both with the BLR visible 
in near infrared. This suggests extended obscuration from the host galaxy is also a plausible 
cause in the absence of PBL in Mrk 334. 

We also note that while using the "T" ratio as indicator of nuclear obscuration, the 
difference between two types of Sy2s becomes more significant (see Fig. 2 (d)). We point 
out that this is mainly because for Compton thick sources, we can only give lower limits of 
~ 10^^ cm"^ for Nh, but the "T" ratio is a continuous variable as long as they are detected 
in the X-ray band. Fig. 7 shows the observed 2-10 keV X-ray flux against the extinction- 
corrected [O III] flux for Compton-thick Sy2s (with a lower limit of 10^"^ cm^^ of Nh, 11 
PBL and 10 NPBL Sy2s). Interestingly, we find although the lower limits of Nh arc the 
same for two types of Compton thick Sy2s, NPBL Sy2s in the figure tend to be weaker in 
X-ray (with a confidence level of 97.6%). We note that two types of Compton-thick Sy2s 
have similar large FeK line EW (> 1 keV) suggesting the X-ray spectra in both types are 
reflection dominated. In this case, smaller "T" ratio in NPBL Sy2s can also be explained by 
higher inclination: sources viewed at higher inclination could have a large fraction of inner 
surface of the torus, where the reflection component is produced, blocked from our line of 
sight, thus expect weaker X-ray emission. 

This work was supported by Chinese NSF through NSF10473009/NSF10533050, and 

the CAS "Bai Ren" project at University of Science and Technology of China. The work 
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the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with NASA. This study has 
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A. Notes on Individual Objects 

In this section, we present brief discussions for the X-ray data of 3 sources in the sample. 
All the errors quoted are at the 90% level of confidence. 

NGC 34: we model the XMM — Newton 0.5 - 8.0 keV spectrum of this source with 
pexrav model (C/dof = 362/448) in terms of its low "T" ratio of 0.03 (<0.1). However, the 
Fe K line is marginally detected with a upper limit EW of 321 eV. We then fit the spectrum 
with an absorbed powerlaw model which gives T = 1.68^q;}7 but no intrinsic absorption, plus 
Fe K fine with EW = 386 (< 1047) eV. The model yields a worse fit with C/dof = 389/448. 
We note the steep spectrum slope may attribute to the host galaxy thermal emission in the 
soft X-ray band. From the lower "T" ratio of 0.03 (<0.1) and better C statistic of pexrav 
model, we consider the galaxy as Compton-thick and give a lower limit of 10^^ cm ~^ to Nh. 

NGC 5728: The Chandra 0.5 - 8 keV spectrum is parameterized here with an absorbed 
powerlaw (F = 2.73±0.17, Nh = T.StH x 10^^ cm^) plus a 0.4% scattered component. The 
Fe K fine is detected at 6.4 keV with EW = 1100^^?^ eV. For the lower "T" ratio and 
large Fe line EW, we then fit the spectrum of this galaxy with pexrav model plus Gaussian 
line. However, the fitting is unacceptable (C/dof = 710/300). So we do not regard it as a 
Compton-thick one in this paper. We note considering it as Compton-thick does not affect 
our results presented here. 

NGC 6552: We fit the XMM - Newton 0.5 - 8 keV spectrum of this source by an 
absorbed powerlaw (F = 2.81q'^3) with a 0.75% scattered component. The best fit (C/dof = 
121/166) gives Nh = 7.i;i:o x 10^^ cm^. The Fe K line is detected at 6.4 keV with EW = 
1408^333 eV. However, the pexrav model can also give an acceptable fitting of the spectrum 
with F = 2.86^[ji^ and Fe K line EW = 4990^1^^ (C/dof = 128/166). The fitted 2-10 keV 
flux is 2.32x10"^^ erg s~^ cm~^. We adopt the absorbed powerlaw model for the spectrum 
fltting in the paper in terms of the better C statistic. We note that the consideration of it 
as Compton thick will not affect our results presented in the paper. 
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Fig. 1. — Rcdshift vs. [O ill] A5007 luminosity for two types of Sy2s in the sample, solid 
circles stand for PBL Sy2s and open circles for NPBL Sy2s. 
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Fig. 2. — The left bottom panel plots the absorption column density Nh versus the [O III] 
A5007 luminosity for PBL and NPBL Sy2s. The distribution of the [O ill] A5007 luminosity 
for two types of Sy2s is presented in the left top panel (solid line for PBL Sy2s and dashed line 
for NPBL Sy2s). The right panel is the plot of FeK line EW (top) and "T" ratio (bottom) 
versus the [O ill] A5007 luminosity. . 
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Fig. 3. — Nh distributions for all the Sy2s in our sample (left) and luminous sources with 
log L[o III] > 41 erg s~^ (right). Shaded areas denote lower limits. The arrows denote the 
upper limits of Nh. 
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Fig. 4. — Distributions of the "T" ratio for all Sy2s (left) and for luminous Sy2s only (right) 
in our sample. Shaded areas denote the 8 sources with only hard X-ray upper limits. 
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Fig. 5. — The plot of "T" ratio vs. extinction-corrected [O ill] flux. Sources with luminosity 
L[o III] <10'^^ erg s~^ are plotted as triangles. 
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Fig. 6. — The plot of "T" ratio vs. [O ill] A5007 luminosity for different spectropolarimetric 
subsamples. (a): Lumsden et al. (2001) 's sample, (b): Tran (2001) 's sample, (c): Moran et 
al. (2000) 's sample, (d): Young et al. (1996) 's sample, (e): other surveys. 
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Fig. 7. — The observed [O ill] flux (extinction corrected) against X-ray (2 - 10 kcV) flux 
for Compton-thick Sy2s (with Nh > 10^^ cm~^). The dashed hnes represent F[o m] = 10 
F2-iofeey (upper) and F[o m] = F2-iofcey (lower). 
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48 


0.69 


>24.126 


2200+fo° 


8 


NGC 3081 


0.0079 


y 


23K 


1.95 


41.43 


27 


1.3 


23.819 


61011?° 


49 


NGC 3281 


0.0115 


n 


23 


1.0 


41.47 


64 


2.9 


24.197 


11801^3°? 


50 


NGC 3982 


0.0037 


n 


15L,23 


0.66 


40.3 


15,31 


0.057 


> 24 


6310+3600 


13 


NGC 4117 


0.0031 


n 


23 








<23.2 






39 


NGC 4388 


0.0084 


y 


lA 


4.51 


41.85 


22 


7.62 


23.43 




5 


NGC 4507 


0.0118 


y 


23K 


4.98 


42.19 


43 


12.8 


23.643 


117 


51 


NGC 4941 


0.0037 


n 


23 


4.57 


41.14 


38 


0.66 


23.653 


16001^°° 


40 


NGC 5135 


0.0137 


n 


2,33A 


6.61 


42.44 


27 


0.16 


>23.954 


1700l«°° 


8 



Table 1 — Continued 



Name 


z 


PBL? 


References 


Fa5007 


L[o III] 


References 


F2-IO keV 


logio Nh 


EW(Fe) 


References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 5283 


0.0104 


n 


15L,23 


0.4 


40.98 


22 


1.46 


23.176 


<220 


7 


NGC 5728 


0.0094 


n 


23, lA 


7.61 


42.18 


27 


1.33 


23.89 


1100+o?n 


13 


NGC 5643 


0.004 


n 


23 


6.62 


41.37 


43 


0.84 


23.845 


500 


52 


NGC 5347* 


0.0078 


y 


42K 


1.14 


41.19 


15 


0.22 


>24 


1300±500 


63 


NGC 5695 


0.014 


n 


15L,23 


0.081 


40.55 


22 


<0.01 






39 


NGC 5929* 


0.0083 


y 


42K 


1.53 


41.40 


2 


1.35 


22.629 




9 


NGC 6890 


0.0081 


n 


23 


0.5 


40.86 


43 










NGC 7672 


0.0134 


n 


20L 








28.6 






39 


Young et al. (1996) 's sample: 






















IRAS 00521-7054 


0.0689 


n 


lA 


0.36 


42.62 


1 


<31.8 






39 


IRAS 04103-2838 


0.118 


n 


lA 
















IRAS 04210+0400 


0.046 


n 


lA 


0.554 


42.44 


1 










IRAS 04229-2528 


0.044 


n 


lA 


0.216 


41.99 


1 










IRAS 04385-0828 


0.0151 


y 


15LP 


0.086 


40.64 


32 


2.4 






39 


IRAS 05189-2524 


0.0426 


y 


lA 


1.3 


42.74 


2 


4.3 


22.756 


30t|g 


2 


IRAS 11058-1131 


0.055 


y 


lA 


0.394 


42.45 


1 


0.39 


>24 


900 


18 


MCG -3-34-64 


0.0165 


y 


lA 


4.0 


42.39 


2 


4.0 


23.614 


35611^3 


67 


IRAS 08277-0242 


0.041 


n 


lA 


1.42 


42.75 


16 








IRAS 13452-4155 


0.039 


n 


lA 
















ESQ 273-IG04 


0.039 


y 


lA 


0.85 


42.48 


1 










IRAS 15480-0344 


0.03 


y 


15P,1A 


5.03 


43.02 


16 


0.37 


>24.204 


<2400 


7 


IRAS 20210+1121 


0.056 


n 


lA 


2.73 


43.31 


1 


0.29 


>25 


1650 


18 


IRAS 20460-1-1925 


0.181 


y 


lA 


0.112 


43.02 


1 


1.5 


22.398 


260+J^I 


49 


IRAS 22017+0319 


0.0611 


y 


15P,1A 


0.42 


42.58 


1 


3.6 


22.69 


380ll«°o 


18 


IRAS 23128-5919 


0.045 


n 


lA 


0.101 


41.68 


1 


0.13 


22.681 




53 


NGC 5506 


0.0062 


n 


15 


3.33 


41.45 


58 


58 


22.46 




11 


Mrk 463E 


0.051 


y 


20L,1A 


1.25 


42.89 


22 


1.46 


23.51 


340+1° 


6 


NGC 4388 


0.0084 


y 


lA 


4.51 


41.85 


22 


7.62 


23.43 




5 


NGC 5252 


0.023 


y 


lA 


0.921 


42.05 


1 


10.7 


22.461 


44±28 


13 


NGC 5728 


0.0094 


n 


23, lA 


7.61 


42.18 


27 


1.33 


23.89 


1100t|?° 


13 


NGC 7496 


0.005 


n 


lA 


0.3 


40.22 


1 


<8 


22.699 




39 


NGC 7674 


0.0289 


y 


1A,20L 


1.93 


42.57 


22 


0.7 


> 24 


0-0+160 
"^'^-170 


14 



Table 1 — Continued 





2 


PBL? 


References 




-r 

^[O III] 


R,efer6nc6S 


^ 2. — iu keV 


ln0"i n IVu 


EW(Fe) References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


Other surveys: 




















Mrk 1210 


0.0135 


y 


54L 


5.8 


42.37 


55 


9.3 


23.263 


108lgg 56 


IRAS 18325-5926 


0.0202 


y 


58A 


1.68 


42.19 


58 


10 


22.31 


242 30 


MCG -5-23-16 


0.008 


y 


58A 


4.09 


41.81 


49 


70 


22.25 


35.2tfg^ 59 


Circinus 


0.0014 


y 


29 


19.1 


40.92 


17 


14 


24.633 


2250l2f» 4 


Mrk 477 


0.038 


y 


54L 


12.4 


43.62 


22 


1.2 


>24 


490^250 49 


NGC 2992 


0.0077 


y 


58A 


1.49 


41.3 


58 


4.5 


21.84 


514±190 49 


NGC 7314 


0.0047 


y 


58A 


17.7 


42.41 


49 


41.2 


22.02 


14711^8 60 


NGC 6300 


0.0037 


n 


58A 


3.2 


40.99 


58 


21.6 


23.342 


148tjg 65 


NGC 7212 


0.027 


y 


54L 


3.2 


42.73 


43 


0.69 


>24.204 




NGC 7590 


0.0053 


n 


33A 


0.168 


40.02 


27 


1.2 


<20.964 


49 


Was 49b 


0.063 


y 


54L 


33.8 


42.51 


57 


0.63 


22.799 


620 ± 250 21 



Note. — The Sy2s with footnote * denote the PBL was detected in later spectropolarimctric observation. Telescope: C = CTIO (4m), P = 
Palomar (5m), K = Keck (10m), L = Lick (3m), S = Subaru (8.2m), E = ESQ (3.6m), A = AAT (3.9m). 

References. — (1) Young et al. 1996; (2) Lumsdcn et al. 2001; (3) Antonucci ct al. 1985; (4) Smith & Wilson 2001; (5) Cappi et al. 2006; (6) 
Imanishi & Terashima et al. 2004; (7) Guainazzi et al. 2005a; (8) Guainazzi et al. 2005b; (9) tartarus.gsfc.nasa.gov; (10) Inglis et al. 1993; (11) 
Bianchi et al. 2003; (12) Levenson et al. 2005; (13) This work; (14) Bianchi et al. 2005a; (15) Tran 2001; (16) de Grijp et al. 1992; (17) Oliva et 
al. 1994; (18) Ueno et al. 2000; (19) Colina et al. 1991; (20) Miller et al. 1990; (21) Awaki et al. 2000; (22) Dahari & De Robertis 1988; (23) 
Moran et al. 2000; (24) Murayama et al. 1998; (25) Matt et al. 2003; (26) Tran 1995; (27) Storchi-Bergmann et al. 1995; (28) Cruz-Gonzalez 
et al. 1994; (29) Alexander et al. 2000; (30) Iwasawa ct al. 2004; (31) Ho ct al. 1997; (32) Tran 2003; (33) Heisler ct al. 1997; (34) Veillcux 
et al. 1995; (35) Vaceli et al. 1997; (36) Prieto ct al. 2002; (37) de Robertis & Ostcrbrock 1986; (38) Storchi-Bcrgmann & Pastoriza 1989; (39) 
Pollctta et al. 1996; (40) MaioUno ct al. 1998; (41) Whittle et al. 1992; (42) Moran et al. 2001; (43) Mulchacy ct al. 1994; (44) Acker ct al. 1991; 
(45) Bianchi et al. 2005b; (46) Kay & Moran 1998; (47) PraqueUi et al. 2003; (48) Lonsdale et al. 1992; (49) Bassani et al. 1999; (50) Vignali & 
Comastri 2002; (51) Matt et al. 2004; (52) Guainazzi et al. 2004; (53) Pranceschini et al. 2003; (54) Tran et al. 1992; (55) Terlevich et al. 1991; 
(56) Masanori et al. 2004; (57) Moran et al. 1992; (58) Lumsden et al. 2004; (59) Dewangan et al. 2003; (60) Dewangan & Griffiths 2005; (61) 
KoUatschny et al. 1983; (62) Risaliti et al. 2000; (63) Levenson et al. 2006; (64) Storchi-Bergmann et al. 1992; (65) Matsumoto et al. 2004; (66) 
Ruiz et al. 1994; (67) Dadina & Cappi 2004; (68) Due et al. 1997; (69) Pernechele et al. 2003; (70) Braito et al. 2003. 



Table 2. Best-fit parameters for the X-ray spectral analysis 



Name 


XMM/Chandra obs. date 


Nh 


r Center Energy 


EW(Fe K) 


/,(%) C/dof F2_l0fceV' 


T 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) (8) (9) 


(10) 



NGC 34 


XMM 2002 Dec 22 


>100 


2.38±0.1 


NGC 3982 


XMM 2004 Jun 15 


>100 


3.74lJ;fj 


NGC 5728 


Chandra 2003 Jun 27 


78i;i 


2.73±0.17 


NGC 6552 


XMM 2002 Oct 18 


nits 


o Q+0.37 
^•^-0.13 


NGC 7172 


XMM 2002 Nov 18 


8.7±0.57 


1.49±0.13 


Mrk 573 


XMM 2004 Jan 15 


>100 


O y + O.OS 

' -0.10 


IRAS 01475-0740 


XMM 2004 Jan 21 


'-'■■-'^-0.02 


2.06±0.06 


Mrk 1066 


Chandra 2004 Jul 18 


>100 


9 7C+0.17 
^- '"^-0.07 



6.4t 


<321 




362/448 


0.23 


0.03 


6.4t 


(,010+3500 




109/135 


0.057 


0.086 


6.39±0.03 


11001320 


0.4±0.03 


299/300 


1.33 


0.17 


6.4t 


1408l««« 


'''-0.22 


121/166 


0.43 


0.27 


6.4t 


40±30 


0.54±0.1 


463/495 


22 


550 


6.4t 


28001 1«20 




394/399 


0.12 


0.07 


6.4t 


130(<344) 




631/776 


0.75 


0.92 


6 34+0- 1'^ 
D-J'*_0.06 


iiZU_ggg 




497/505 


0.23 


0.05 



Note. — Col. (1): galaxy name. Col. (2): Observation date. Col. (3) power-law photon index. Col. (4); measured absorption column density, in units 
of 10^2 cm~2. Col. (5): the Fe line energy in units of keV. Col. (6): the Fe line equivalent width in units of eV. Col. (7): scattering fraction of the soft 
component. Col. (8): C statistic and number of degrees of freedom (dof). Col (9): the fitted 2-10 Kev flux in units of 10"^^ erg cm"^. Col (10) T: the 
ratio of F2_io kev/F[0 iii]- ^ : fixed. 
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Table 3. A summary of the Statistical Properties for Sy2s 



Parameters 


Sy2s with PBL 


Sy2s without PBL 


Note 


PSuII 


log(T) 


-0.087±0.145 


-0.342±0.217 


total 


25% 


log{T) 


-0.201±0.138 


-0.87±0.186 


luminous 


0.9% 


logCT)*- 


-0.093±0.144 


-0.436±0.199 


total 


13.5% 


logCT)*- 


-0.203± 0.137 


-0.899±0.179 


luminous 


0.68% 


log[EW{Fe)] 


2.654±0.107 


2.85±0.127 


total 


18.6% 


log[EW(Fe)] 


2.626±0.107 


2.999±0.066 


luminous 


4.67% 


log(NH) 


23.755±0.190 


23.852±0.274 


total 


66.5% 


log(NH) 


23.739±0.212 


24.428±0.192 


luminous 


7.7% 


log(L[o III]) 


42.099±0.13 


41.338±0.188 


total 


0.34% 


log(L[o III]) 


42.198±0.121 


41.991±0.157 


luminous 


30.6% 


log(z) 


0.020±0.003 


0.017±0.004 


luminous 


29.4% 



Note. — " the possibility Pnull is for the null hypothesis that the two dis- 
tributions are drawn at random from the same parent population. ^ including 

8 hard X-ray upper limits stated in §2. When there are censored data, wc use 
Gehan's generalized Wilconxon test-permutation variance (GGW test, one kind 
of ASURV test), "total": all Sy2s with X-ray data in the sample, "luminous": 
for Sy2s with Ljq m] > 10'*^ erg only. 



